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For  many  types  of  superconducting  qubits,  magnetic  flux  noise  is  a  source  of  pure  dephasing.  Measurements 
on  a  representative  dc  superconducting  quantum  interference  device  (SQUID )  over  a  range  of  temperatures  show 
that  S$(f)  =  A2/(//l  Hz)“,  where  5*  is  the  flux  noise  spectral  density,  A  is  of  the  order  of  1  p.<t>o  Hz-1/2, 

0.61  ^  a  ^  0.95,  and  (ho  is  the  flux  quantum.  For  a  qubit  with  an  energy  level  splitting  linearly  coupled  to  the 
applied  flux,  calculations  of  the  dependence  of  the  pure  dephasing  time  of  Ramsey  and  echo  pulse  sequences 
on  a  for  fixed  A  show  that  decreases  rapidly  as  a  is  reduced.  We  find  that  is  relatively  insensitive  to  the 
noise  bandwidth,  /i  ^  ^  f2,  for  all  a  provided  the  ultraviolet  cutoff  frequency  /2  >  l/r^,.  We  calculate  the 

ratio  of  the  echo  (£)  and  Ramsey  (R)  sequences  and  the  dependence  of  the  decay  function  on  o'  and 

f2.  We  investigate  the  case  in  which  S<j>(/o)  is  fixed  at  the  “pivot  frequency”  /0  ^  1  Hz  while  a  is  varied  and  find 
that  the  choice  of  fo  can  greatly  influence  the  sensitivity  of  and  r ^  to  the  value  of  a.  Finally,  we  present 
calculated  values  of  in  a  qubit  corresponding  to  the  values  of  A  and  a  measured  in  our  SQUID. 

DOI:  10.1 103/PhysRevB. 85.224505  PACS  number(s):  85.25.Dq,  03.67.Lx,  05.40.Ca 


I.  INTRODUCTION 

The  dynamics  of  superconducting  quantum  bits  (qubits)1 — 
broadly  classified  as  charge  qubits,2  flux  qubits,’  and  phase 
qubits4 — can  be  characterized  by  two  times:  the  relaxation 
time  T\  and  the  pure  dephasing  time  r^.5  The  time  7) 
required  for  a  qubit  to  relax  from  its  first  excited  state  to  its 
ground  state  is  determined  by  the  strength  of  environmental 
fluctuations  at  a  frequency  corresponding  to  the  energy  level 
splitting  t>oi  of  the  two  states.  The  decoherence  time  73,  over 
which  the  phase  of  superpositions  of  two  eigenstates  becomes 
randomized,  has  two  contributions:  1/72  =  1/(2  7))+  l/r^. 
The  pure  dephasing  time  r,,,  is  limited  by  fluctuations  in  voi, 
due  predominantly  to  fluctuations  in  magnetic  flux  in  the  case 
of  flux  qubits. 

Excess  low-frequency  flux  noise  was  first  identified  in  dc 
superconducting  quantum  interference  devices  (SQUIDs).6 
Measurements  at  millikelvin  temperatures7-10  reveal  a  power 
spectrum  Sq,(f)  scaling  as  1  //“  (/  is  frequency),  with  an 
amplitude  at  1  Hz  typically  of  the  order  of  1  fi<t> o  Hz-1/2,  that 
is  surprisingly  uniform  for  SQUID  washers  of  widely  differing 
geometries  that  are  fabricated  with  a  variety  of  materials.  Here 
cj>0  =  h /2e  is  the  flux  quantum. 

Flux  noise  is  believed  to  arise  from  the  random  reversal 
of  electron  spins  at  the  interface  between  a  superconducting 
film  and  an  insulator.11-13  The  areal  density  of  independent 
spins  required  to  account  for  the  observed  flux  noise  is 
about  5  x  1017m-2,  a  value  that  has  been  corroborated  by 
observations  of  paramagnetism  in  SQUIDs9  and  normal  metal 
rings.14  Recent  experiments  on  anticorrelations  of  flux  noise 
have  confirmed  the  surface  spin  model.1516  An  unambiguous 
understanding  of  the  mechanism  by  which  the  spins  produce 
1  /f  flux  noise,  however,  has  yet  to  be  developed.  Recent 


proposals  include  spin  clusters,17  spin  glasses,18  fractal  spin 
clusters,19  and  hyperfine  interactions;20  the  models  in  Refs.  18 
and  1 9  suggest  that  a  may  differ  from  unity. 

Measurements  of  xv  in  flux  qubits15,16,21-23  and  phase 
qubits24  have  been  used  to  infer  the  magnitude  of  the  flux 
noise  in  these  devices,  under  the  assumption  that  the  spectral 
density  of  the  flux  noise  scaled  as  1//.  In  this  paper  we  first 
present  measurements  of  flux  noise  spectral  densities  scaling 
as  1  //“  in  which  the  exponent  a  deviates  markedly  from 
unity.  We  then  show  theoretically  that  such  deviations  strongly 
affect  tv:  The  value  of  T,f  decreases  markedly  with  decreasing 
a.  Additionally,  we  examine  the  influence  of  a  and  noise 
bandwidth  on  r^,  the  ratio  r^,  r/t^.r  obtained  in  echo  ( E )  and 
Ramsey  ( R )  pulse  sequences,  and  the  functional  dependence 
of  the  decay  function.  Finally,  we  calculate  the  predicted 
for  values  of  A  and  a  obtained  in  our  measurements. 


II.  EXPERIMENTAL  PROCEDURES  AND  RESULTS 

We  measured  the  flux  noise  spectral  densities  of  Nb-based 
dc  SQUIDs,  fabricated  using  a  50-A/cm2  Nb/A10.v/Nb  trilayer 
process.25  Each  2.5  x  2.5-//m2  junction  was  shunted  with  a 
2.5-G  PdAu  resistor  to  eliminate  hysteresis  on  the  current- 
voltage  (I-V)  characteristic.  Up  to  six  SQUIDs,  connected  in 
series,  were,  in  turn,  connected  in  series  with  a  compensating 
resistor  Rc  ~  0.45  G  and  the  superconducting  input  coil  of  a 
readout  SQUID,  operated  in  a  flux-locked  loop  (Fig.  1).  To 
measure  the  noise  in  a  given  SQUID,  we  applied  a  current 
lb  sufficient  to  produce  a  voltage  V  of,  typically,  5  //V 
across  it.  The  resulting  static  current  around  the  circuit  was 
canceled  by  an  appropriate  current  /,  in  R,  to  ensure  that 
the  remaining  SQUIDs  remained  in  the  zero-voltage  state.  In 
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FIG.  1.  (Color  online)  Circuit  schematic  of  measurement  system 
for  three  measured  SQUIDs.  Measurement  of  the  middle  SQUID  is 
shown  (see  text).  A  single  coil  applies  flux  to  all  SQUIDs. 


addition,  a  choke  inductor  was  used  to  decouple  oscillations 
at  the  Josephson  frequency  between  the  measured  and  readout 
SQUIDs.  Since  Rc  was  much  less  than  the  dynamic  resistance 
of  any  given  SQUID,  the  SQUID  was  effectively  voltage 
biased.  Fluctuations  in  the  critical  current  of  the  measured 
SQUID  induced  a  current  noise  I  with  spectral  density  Sj(f ) 
in  the  input  coil  of  the  readout  SQUID.  We  inferred  the 
flux  noise  from  =  S/(/)/[(3//3<I>)y]2,  where  is 

the  flux  applied  via  an  external  coil  to  the  measured  SQUID 
and  (3//3<J))y  was  determined  separately.  We  also  measured 
the  critical  current  noise  of  the  junctions6  in  each  SQUID 
biased  at  («  is  an  integer)  so  (3//3<t>)y  =  0.  This  noise 
was  negligible  compared  with  the  current  noise  produced  by 
the  flux  noise  for  large  values  of  ( 3  /  /  3  <t> )  y .  The  experiment 
was  mounted  in  a  lead-coated  copper  box  surrounded  by  a 
cylindrical  lead  shield  inside  a  cryoperm  shield  and  cooled 
with  a  dilution  refrigerator.  All  leads  were  heavily  filtered. 

Figure  2  shows  power  spectra  of  a  single  SQUID  at  three 
different  temperatures.  The  inner  and  outer  dimensions  of 
the  washer  were  50  and  90  gm,  respectively.  At  higher 
frequencies,  the  spectra  begin  to  flatten  out  due  to  white  current 
noise  from  the  SQUID  shunt  resistors,  which  dominates  that 
from  Rc.  The  spectra  were  fitted  to  the  form 

S(f)  =  A2/(//l  Hz)"  +  C2  (1) 


/(Hz) 

FIG.  2.  (Color  online)  Measured  and  fitted  flux  noise  spectra. 
Measured  spectra  at  0.2,  1.3,  and  4.5  K  are  continuous  curves; 
symbols  are  fits  to  Eq.  (1).  Quoted  values  in  the  fits  are  in  units 
of(/z<I>o)2/Hz. 


with  parameters  A  and  C  for  the  amplitude  of  the  “1//”  flux 
noise  and  white  noise,  respectively.  We  found  that  the  exponent 
a  of  the  1  //"  noise  can  be  far  from  unity,  varying  from  0.61 
to  0.95.  We  remark  that  these  values  are  representative  of 
measurements  on  about  20  SQUIDs.  Since  the  junctions  in  flux 
qubits  are  not  resistively  shunted,  we  shall  focus  on  dephasing 
from  the  term  A2/(//l  Hz)". 


III.  THEORETICAL  RESULTS 


A.  Model 

Since  a  can  evidently  be  much  less  than  unity,  it  is 
natural  to  ask  what  impact  this  has  on  the  pure  dephasing 
of  flux  qubits.  Low-frequency  flux  noise  modulates  the  energy 
splitting  of  the  ground  and  first  excited  states  of  a  flux  qubit, 
hv oi  =  [A2  +  e2(d>)]1/2,  via  the  bias  energy  e(<J>).  The  bias 
energy  is  the  energy  difference  between  the  two  states  with 
persistent  currents  ±/?  when  there  is  no  tunneling  between 
them  (A  =  0).26  Here  e  —  2 Iq ( <t>  —  <t>0/2)  or,  equivalently, 
1 q  =  |(3e/3d>). 

We  define  the  sensitivity  of  the  splitting  to  a  change  in  <f>  in 
terms  of  the  longitudinal  sensitivity  of  the  qubit  to  flux  noise, 

=  3uoi/3<F  =  (l//?)(3e/3<I>)e/(A2  +  e2)1/2.  (2) 


To  first  order,  there  is  no  dephasing  from  flux  noise  at  the 
degeneracy  point,  <t>  =  <t>o/2,  where  e  vanishes.  In  this  paper, 
however,  we  consider  the  limit  e/A  1,  far  from  the  de¬ 
generacy  point,  at  which  Svoi/SO  =  (l//i)3e/3<£>  —  2 Iq/ h. 
We  assume  that,  in  this  limit,  l/r^  l/27j  so  the  measured 
dephasing  time  arises  only  from  pure  dephasing.  We  adopt 
the  value  D®  =  1012  Hz/<t>o.  corresponding  to  the  typical 
value21,22  Iq  ~  0.3  //A.  Furthermore,  based  on  the  empirical 
observation  that  S<j>(1  Hz)  is  relatively  constant  among  a 
wide  variety  of  SQUIDs,  we  assume  that  A  —  1  //To  Hz”1/2 
regardless  of  the  value  of  a.  We  consider  noise  fixed  at  a 
frequency  other  than  1  Hz  in  Sec.  Ill  F. 

The  modulation  of  v>oi  by  flux  noise  leads  to  an  accu¬ 
mulation  of  phase  error  and,  thus,  to  dephasing.  The  rate  at 
which  the  dephasing  occurs  varies  between  different  types  of 
pulse  sequences.  For  example,  in  a  Ramsey  sequence27  the 
qubit  is  excited  by  a  microwave  jt/2  pulse  from  the  ground 
state  into  a  superposition  of  ground  and  excited  states.  After 
a  time  t  another  7r/2  pulse  is  applied  and  the  qubit  state  is 
measured.  The  results  of  many  measurements  with  fixed  t  are 
averaged  and  t  is  varied  from  t  T,p  to  t  'A>  r,/;  to  obtain  the 
decay  function  g(t).  Here,  we  define  the  dephasing  time  as 
gitfj,)  =  1/e.  To  eliminate  dephasing  due  to  flux  fluctuations 
between  pulse  sequences,  one  implements  an  echo  sequence 
in  which  a  tt  pulse  is  inserted  midway  between  the  two  jt/2 
pulses.28  In  general,  the  echo  sequence  yields  a  dephasing  time 
greater  than  that  of  the  Ramsey:  ^  >  t(/kK. 

The  sensitivity  of  the  Ramsey  and  echo  sequences  to  noise 
are  described  by  the  weighting  function  W(f,t )  given  by4'29 


sin2(7r/r) 

(tt  ft)2 


WE(f,t ) 


sin4(7r  ft/ 2) 

(nft/ 2)2 


(3) 


For  the  Ramsey  sequence  with  /  <(C  1/f,  we  see  that  WR  m  1, 
whereas  for  /  1/f,  Wr  falls  as  l//2.  Consequently,  we 

expect  the  dominant  contributions  to  the  Ramsey  dephasing 
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FIG.  3.  (Color  online)  Computed  values  of  r and  t0>e  vs.  a  for 
/i  =  10-3,  0.1,  10,  and  103  Hz  and  f2  oo. 


time  to  arise  from  noise  at  frequencies  /  <  I  /t0iR.  In  contrast, 
for  the  echo  sequence  WE  scales  as  f2  for  /  \/t  and  as  1  // 2 

for  /  »  1/L  In  this  case,  we  expect  the  dominant  contribution 
to  the  dephasing  time  to  be  from  noise  at  frequencies  /  ~ 
1  /?>,£. 

The  decay  function  g(t)  is  calculated  by  ensemble  averag¬ 
ing  over  the  entire  measurement  time,  yielding4-29 


g(t)  =  exp 


—t2(2n  D<t,)2 


(4) 


Here,  the  symmetrized  noise  power  is  defined  as  S®(f)  = 
(l/2)/off{(O(r)O(0))  +  (tp(0)d>(f))}e_2lr'^,  which  we  re¬ 
place  with  the  observed  spectrum:  S®(/ )  =  A2/[//(l  Hz)]"; 
/i  and  /2  are  cutoff  frequencies  limiting  the  noise  frequency 
bandwidth  to  which  the  qubit  is  sensitive.  Independent  of  the 
particular  pulse  sequence,  the  infrared  cutoff  f\  is  set  by 
the  entire  measurement  time  T  taken  to  acquire  sufficient 
statistics  to  determine  the  decay  function  g(t),  that  is,  f\  — 
l/T,  where  T  may  range  from,  say,  1  ms  to  1000  s.  What 
determines  the  ultraviolet  cutoff  /2,  however,  is  less  clear. 
Recent  experiments30'31  indicate  that  flux  noise  can  not  only 
extend  to  very  high  frequencies  (in  one  case  in  excess  of  1  GHz) 
but  also  maintain  its  nonunity  value  of  a  out  to  /2. 


B.  Dephasing  times  versus  a 

As  is  evident  from  Eq.  (4),  a  nonunity  value  of  a  will  affect 
the  integral  in  a  complicated  way.  Figure  3  shows  computed 
dephasing  times  for  both  sequences  versus  a  for  /2  ->  oo  and 
f\  —  1  (L3,  1  0“  1 ,  101,  and  103  Hz.  The  effect  of  changing 
a  is  substantial:  both  r0  R  and  x,pE  increase  by  an  order  of 
magnitude  as  a  is  varied  from  0.6  to  0.9.  By  comparison, 
we  find  that  an  order  of  magnitude  change  in  A  for  a  given 
value  of  a.  also  changes  r0  by  an  order  of  magnitude.  Figure  3 
further  shows  that,  because  of  its  insensitivity  to  low-frequency 
noise,  the  echo  sequence  yields  significantly  longer  dephasing 
times  for  all  a.  Finally,  while  r0i£  is  insensitive  to  changes 
in  /i  for  /i  «;  l/r0,£  (equivalently  T  »  r0j£),  r 0,R  becomes 
increasingly  sensitive  as  a  increases. 


C.  Dephasing  times  versus  cutoff  frequencies 

We  now  examine  more  quantitatively  the  sensitivity  of  r0 
to  changes  in  both  f\  and  /2  for  various  values  of  a.  For 


FIG.  4.  (Color  online)  Normalized  Ramsey  dephasing  times  for 
0.6  ^  a  ^  1.2  in  steps  of  0.1.  (a)  r0 ,r(/i)/t0iR(/i  =  0.1  Hz)  versus 
/i  for  f2  -*■  oo  and  (b)  t0iR(/2)/t0iR(/2  -*■  oo)  versus  f2.  The  dots 
in  (b)  are  placed  at  f2  =  1  /r0,R(/2  -*■  oo),  above  which  displays 
no  dependence  on  f2  (see  text). 

the  Ramsey  sequence  with  /2  -»  oo.  Fig.  4(a)  shows  r ^  R, 
normalized  to  =0.1  Hz)  versus  f\  for  0.6  ^  a  ^  1.2. 

We  again  see  that  the  sensitivity  of  ^  to  f\  increases  with 
increasing  a.  Even  so,  for  a  —  1.2,  r  changes  by  a  factor 
of  only  4  when  f)  is  varied  from  0.1  to  104  Hz. 

To  explore  the  effect  of  f2  on  we  fix  f\  =  1  Hz  and 
vary  /2,  plotting  r0^(/2)/r0,R(/2  ->  oo)  for  0.6  ^  a  <  1.2 
[Fig.  4(b)].  We  see  that  the  sensitivity  of  r0,R  to  /2  increases 
for  decreasing  a.  Furthermore,  Fig.  4(b)  shows  that  r0  R  is 
insensitive  to  the  particular  value  of  /2  for  /2  1/t0,r(/2  — > 

oo),  simply  because  the  Ramsey  sequence  is  insensitive  to 
noise  for  /  l/r0  R.  However,  as  /2  decreases  through 
l/r0,R(/2  —*■  oo)  a  non-negligible  amount  of  noise  to  which 
the  qubit  is  sensitive  is  effectively  eliminated,  thereby  reducing 
the  total  integrated  noise  and  increasing  r0,R(/2).  This  effect 
is  greater  for  small  a,  where  S®  decreases  with  /  more  slowly 
and  contributes  to  dephasing  out  to  a  higher  frequency. 

We  perform  a  similar  analysis  of  the  sensitivity  of  r0  R  to  the 
value  of  /2.  In  Fig.  5  we  plot  r0,£(/2)/r0i£(/2  oo)  versus 
/2  for  /|  =  I  Hz.  As  with  the  Ramsey  sequence,  we  find  that 
r0i£  is  insensitive  to  /2  for  /2  l/r0i£(/2  -*  oo).  Indeed, 

since  r0iR  is  dominated  by  noise  at  /  ~  1  /r0>R,  this  result  is  as 
we  expect.  Also  in  analogy  with  the  Ramsey  sequence,  r0iR  is 
more  sensitive  to  /2  for  small  a.  Unlike  the  Ramsey  sequence, 
however,  where  the  dephasing  is  sensitive  to  frequencies  over 
a  large  bandwidth  (f  \  to  1  /r0  R),  the  echo  sequence  is  sensitive 
to  noise  only  in  a  narrow  bandwidth  around  l/r0,£,  making 
r0  £•  much  more  sensitive  to  changes  in  /2  for  /2  ^  l/r0  E- 
Here,  r0,R  increases  by  an  order  of  magnitude  for  a  two-order- 
of-magnitude  decrease  in  /2. 

D.  The  ratio  r0,£/r0  R 

Since  the  value  Z)0  can  vary  significantly  between  flux 
qubits,  we  consider  the  ratio  r0i£/r0,R,  which  has  the 
advantage  of  being  rather  insensitive  to  the  precise  values  of 
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h  (Hz) 


FIG.  5.  (Color  online)  Computed  values  of  1 1$ ,E,{fi  — >■ 

00)  vs.  f2.  Lower  cutoff  frequency  ft  =  1  Hz  and  0.6  ^  a  ^  1.2  in 
steps  of  0.1.  The  colored  dots  are  placed  at  f2  =  l/x^^ifi  — >  00), 
above  which  displays  no  dependence  on  f2  (see  text). 

both  A  and  D^.  We  compute  these  times  using  Eq.  (4),  which 
shows  that  the  decay  function  g{t)  depends  only  on  the  product 
A  Dtp.  To  explore  the  dependence  of  the  ratio  on  a,  we  compute 
versus  a  for  f2  00  (equivalent  to  /1  »  1  / r^)  and 
fi  =  10_1,  101,  and  103  Hz.  Furthermore,  for  each  value  of 
f\  we  perform  the  calculation  for  AD^/i  106  Hz1/2)  =  0.2,  1, 
and  5.  The  results  are  shown  in  Fig.  6. 

First,  we  examine  the  dependence  on  a.  As  u  increases, 
noise  at  frequencies  much  greater  than  1  Hz  falls  quickly,  so 
^  increases  rapidly.  Conversely,  noise  at  low  frequencies 
near  1  Hz  changes  little  as  a  changes.  The  Ramsey  dephasing 
time  is  sensitive  to  a  large  noise  bandwidth  where  a  significant 
contribution  comes  from  frequencies  near  f\ .  Therefore,  as  a 
increases  we  expect  r  to  increase  less  rapidly  than  t^e, 
explaining  the  increasing  trend  of  t^e/t^e 

For  small,  fixed  values  of  a,  changing  the  value  of  f\ 
changes  the  ratio  only  slowly  because  both  t^e  and  r,:,  K  are 
limited  by  noise  at  /  ^>  f\ .  As  the  value  of  a  increases,  how¬ 
ever,  an  increasing  contribution  to  dephasing  in  the  Ramsey 
sequence  arises  from  lower  frequencies  /  f\ .  Therefore,  for 
large,  fixed  values  of  a,  increasing  f\  has  the  effect  of  removing 
a  significant  noise  contribution,  thereby  increasing  t^e  and 


a 


FIG.  6.  (Color  online)  Ratio  vs.  a.  Lower  cutoff 

frequency  f\  =  10_1,  101,  and  103  Hz  and  f2  — >  00.  The  thin  upper, 
heavy  middle,  and  thin  lower  lines  correspond  to  AD^/ilO6  Hz1/Z)  = 
0.2,  1,  and  5,  respectively. 


decreasing  the  ratio  t^e/t^e-  We  remark  that  since  f\  is  an 
experimentally  variable  parameter,  measuring  r ^e/t^.e  for 
several  different  measurement  times  may  shed  light  on  the 
value  of  a. 

Finally,  we  see  that  the  r 0,e/t0,e  is  moderately  sensitive 
to  the  product  AD $  only  for  a  >  1.  However,  additional 
calculations  show  that,  for  f2  <  l/r^,  the  ratio  becomes 
extremely  sensitive  to  the  particular  value  of  A  Dp. 

E.  Dependence  of  decay  function  on  of  and 
ultraviolet  cutoff  frequency 

The  decay  function  is  of  particular  interest  experimentally, 
since  it  can  be  measured  directly.  In  general,  the  decay 
function  of  7) -limited  processes  is  a  simple  exponential, 
that  is,  git)  —  exp(— f  /  7] ).  However,  the  decay  function  of 
pure  dephasing  processes  is  more  complicated  and  can  be 
characterized  as  g(t)  =  exp[—  x(0L  where  x(0  can  contain 
terms  that  are  higher  order  in  t. 

Here,  we  examine  the  functional  dependence  of  /(f)  for 
both  pulse  sequences.  In  each  case,  we  find  that  x(t)  oc  tY , 
where  y  can  take  two  values  (y\  and  y2)  within  a  single 
sequence,  separated  by  a  characteristic  time  set  by  \/f2: 
X(t  I//2)  oc  tn  and  x(?  I//2)  oc  tYl.  For  the  Ramsey 

sequence,  y \  —2  and  y2  =  1  +  a  for  a  1  and  y2  =  2  for 
a  >  1.  For  the  echo  sequence,  yi  =  4  and  y2  =  1  +  a.  These 
results  reveal  two  experimentally  relevant  insights.  First,  for 
1  1  If 2,  y2  depends  on  a.  Thus,  if  1  /f2,  a  careful  fit  of 
the  experimentally  observed  decay  envelope  may  shed  light  on 
the  value  of  a.  Second,  the  functional  form  of  g(t)  can  reveal 
information  about  f2.  For  example,  if  one  does  not  observe  that 
X  (f)  oc  f4  in  an  echo  experiment,  f2  must  be  as  high  as  1  / t^e, 
establishing  an  important  lower  bound  on  the  bandwidth  of  the 
flux  noise. 

Figure  7  emphasizes  the  above  statements,  showing  git) 
plotted  for  both  sequences  for  a  —  0.6  and  1.2,  and  for  f2 
both  above  and  below  1  jxp(f2  —*■  00),  thereby  showing  both 
y  1  and  y2  dependence.  In  Figs.  7(a)  and  7(b)  we  plot  g(t) 
for  the  Ramsey  sequence  with  a  —  0.6  and  1.2.  We  note  that 


0  12  12 

t/Tj 

FIG.  7.  (Color  online)  Computed  decay  function  g(t)  versus  t/rp 
for  (a)  and  (b)  Ramsey  sequences  and  (c)  and  (d)  echo  sequences 
with  a  =  0.6  and  1.2.  In  the  dashed  red  trace,  f2  1/t p(f2  — >  00); 
in  the  solid  blue  trace  f2  1/t^(/2  — >■  00). 
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FIG.  8.  (Color  online)  Computed  dephasing  times  vs.  a  for 
f\  =  1  Hz  and  /)  — >  oo.  (a)  Ramsey  and  (b)  echo  pulse  sequences 
for  fixed  S<i,(/o)>  where  the  pivot  frequency  f0  =  10"2,  1,  102,  and 
104  Hz. 

the  difference  between  the  functional  dependencies  of  the  two 
traces  in  Fig.  7(a)  is  slight  and  would  be  nearly  impossible 
to  measure  experimentally.  In  Fig.  7(b)  there  is  no  functional 
difference  since  y1  =  y2.  Figures  7(c)  and  7(d)  show  g(t)  for 
the  echo  sequence  for  a  —  0.6  and  1.2.  The  difference  is  more 
dramatic  since  y\  —  4  is  so  large.  In  this  case,  such  a  difference 
might  be  experimentally  observable. 

F.  .S'o(/)  pivoting  about  fn  #  I  Hz  as  a  is  varied 

As  mentioned  previously,  there  is  no  a  priori  reason  to 
hold  5<t>(l  Hz)  fixed  as  a  is  varied;  the  choice  is  based  on  the 
empirical  observation  that  values  of  £$(1  Hz)  are  relatively 
uniform  across  a  wide  variety  of  devices  and  measured  a.  To 
explore  the  sensitivity  of  our  calculations  to  this  assumption, 
we  calculated  the  dephasing  times  for  both  sequences  versus 
a  for  fixed  £<s(/o),  where  /o  =  1 0  2,  1,  102,  and  104  Hz. 
Conceptually,  the  spectra  can  be  imagined  as  pivoting  as  a 
changes  about  a  fixed  spectral  density  S®(fo)  at  frequency 
/o.  In  order  to  normalize  the  magnitude  of  each  set  of  curves 
corresponding  to  a  particular  /o,  we  choose  as  a  convention 
that  .S'o  ( 1  Hz)  =  A2  when  a  =  1,  regardless  of  the  value  of  /o, 
that  is,  S<j,(/ )  =  (A2 / fo)(f / fo)~a .  This  convention  is  based 
loosely  on  empirical  observation;  it  does  not  significantly 
change  the  dependence  of  on  a,  but  merely  sets  the  absolute 
scale. 

The  results  of  these  calculations,  plotted  in  Fig.  8,  show  a 
dramatic  effect,  both  qualitatively  and  quantitatively,  on  the 
dependence  of  t^r  and  x,ike  on  a.  For  the  Ramsey  sequence 
[Fig.  8(a)],  the  general  trend  of  increasing  t^r  is  significantly 
altered  as  /o  increases  and  even  becomes  nonmonotonic 
for  /o  =  104  Hz.  In  addition,  for  small  values  of  a,  t^r 
increases  dramatically  as  /o  increases.  We  note  that,  because 
of  our  normalization  condition,  the  curves  intersect  at  a  —  1 . 
Calculations  for  the  echo  sequence  are  shown  in  Fig.  8(b), 
which  shows  a  similar  dependence  of  t^.r  on  /o.  For  both 
sequences,  the  dependence  of  x(p  on  a  is  minimal  for  /o  = 


TABLE  I.  Computed  x^  R  and  r  with  fi=  1  Hz  and  /2  — >■  oo 
for  flux  qubits  with  Iq  =  0.3  p  A.  (Upper  section)  Values  of  A  and  a 
from  Fig.  2;  (middle  section)  A  reduced  by  factor  of  10,  a  unchanged; 
(lower  section)  A  set  equal  to  1  Hz~1/2,  a  unchanged. 


A(/z<t>o  Hz-‘/2) 

a 

T/,,r  (ns) 

x^,e  (ns) 

1.78 

0.61 

1.2 

2.5 

1.98 

0.79 

5.6 

15.4 

3.35 

0.95 

8.9 

37.8 

0.178 

0.61 

20.9 

43.2 

0.198 

0.79 

73.3 

202.5 

0.335 

0.95 

99.5 

400.8 

1.0 

0.61 

2.4 

5.1 

1.0 

0.79 

11.9 

33.1 

1.0 

0.95 

31.6 

130.5 

104  Hz,  the  highest  computed  /o.  This  dependence  is  easily 
understood  for  the  echo  sequence,  which  is  sensitive  only  to 
noise  at  f  ra  1/t^^e-  As  /o  approaches  1/t^r,  the  effect  of 
a  eventually  becomes  negligible.  In  fact,  if  /o  were  to  exceed 
l/r0,£>  the  trend  in  a  would  actually  reverse.  The  Ramsey 
sequence,  however,  is  sensitive  to  a  larger  noise  bandwidth 
and  has  a  correspondingly  more  complicated  dependence, 
exhibited  by  its  nonmonotonic  behavior  for  large  values  of  /o. 

G.  Tabulated  dephasing  times 

Finally,  we  use  our  theoretical  prediction  of  a  strong 
dependence  of  the  dephasing  times  on  a  to  calculate  r  and 
x,p  e  for  the  experimental  values  of  A  and  a  shown  in  Fig.  2. 
We  assume  f\  =  1  Hz  and  f2  —*■  oo.  The  results  are  shown  in 
the  upper  section  of  Table  I.  We  see  that,  despite  having  the 
largest  value  of  the  flux  noise  magnitude  A,  the  spectrum  with 
the  highest  value  of  a,  0.95,  yields  the  longest  dephasing  times. 
This  result  emphasizes  a  crucial  point:  simply  lowering  the  flux 
noise  magnitude  A  while  keeping  a  constant  may  not  be  the 
most  effective  avenue  to  increasing  r^.  The  middle  section  of 
Table  I  shows  the  effect  on  t^r  and  x,il  K  of  a  10-fold  reduction 
in  A  for  fixed  a.  The  factors  by  which  x,p  R  and  x^  E  increase 
are  comparable  and  decrease  as  a  increases,  from  about  17 
(a  =  0.61)  to  about  11  (a  =  0.95).  The  values  of  t^r  and 
x(i,e  for  A  =  I  /i  <f>(l  Hz-1/2  are  shown  in  the  lower  section  of 
Table  I.  As  expected,  r  and  x^e  increase  dramatically  as  a 
increases  from  0.61  to  0.95. 

IV.  CONCLUDING  REMARKS 

In  conclusion,  we  have  presented  data  showing  that,  in 
general,  flux  noise  scales  as  1  //“,  where  0.6  <  a  <  1.0.  Our 
subsequent  calculations  show  that  the  predicted  dephasing 
times  r^  R  and  of  a  qubit  are  very  sensitive  to  the 
value  of  a.  As  the  value  of  a  increases,  both  t^.r  and 
increase  dramatically — by  an  order  of  magnitude  in  some 
cases.  Since  experimentally  inferred  values  of  5,<d(1  Hz)  from 
qubit  measurements  have  generally  assumed  that  a  —  1,  a 
nonunity  value  of  a  can  introduce  a  significant  error  into  the 
inferred  value  of  A.  Furthermore,  we  have  shown  that  while  the 
lower  cutoff  frequency  f\  (set  by  the  total  measurement  time) 
does  not  significantly  affect  rf/).  the  upper  frequency  cutoff 
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fi  can  significantly  change  in  a  manner  dependent  on  the 
value  of  a,  particularly  for  the  echo  sequence.  Moreover,  we 
have  shown  that  by  examining  the  directly  measurable  ratio 
and  the  dephasing  function  g(t),  experimentalists 
may  have  a  probe  into  the  values  of  a  and  fa.  Finally, 
the  frequency  at  which  the  flux  noise  spectra  pivot  can 
dramatically  affect  the  sensitivity  of  r((!)  to  a. 

Most  importantly,  these  results  demonstrate  that  lowering 
the  flux  noise  amplitude  is  not  the  only  method  of  increasing 
qubit  dephasing  times.  With  a  more  detailed  understanding 
of  what  sets  a  experimentally — for  example,  the  geometry  of 
the  qubit  washer — it  may  be  possible  to  increase  dephasing 


times  substantially  by  raising  the  value  of  a.  Finally,  we  note 
that  with  straightforward  modification  our  formalism  could  be 
used  to  calculate  dephasing  times  from  critical  current  noise 
and  charge  noise  for  the  case  a  ^  1 . 
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